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We discuss the impact of including corrections beyond single gluon 
exchange in light mesons within the nonperturbative framework of Dyson- 
Schwinger equations (DSE) and Bethe-Salpeter equations (BSE). We do 
this by considering unquenching effects in the form of hadronic resonance 
contributions, notably pion exchange, and by the inclusion of the dominant 
gluon self-interactions to the quark-gluon vertex. Thus we make steps 
towards an ab initio description of light mesons by functional methods. 

PACS numbers: ll.10.St, 11.30.Rd, 12.38.Lg, 14.65.Bt 

1. Introduction 

With the absence of quarks and gluons in the physical spectrum, we are 
forced to probe their interaction at low energies by the study of colourless 
composites of these particles. This information is encoded in the Green's 
functions of our QFT, in particular the four quark-scattering matrix and the 
quark-gluon vertex. Faced with the richness of the hadronic spectrum, it is 
not unreasonable to draw a comparison with the complex nonperturbative 
particulars of these Green's functions. In this talk we will consider the 
simplest bound states - the light mesons - as our probes and investigate how 
their properties depend on information present in the quark-gluon vertex. 
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2. Dyson-Schwinger equations 

To solve the DSEs, shown for the inverse quark propagator and quark- 
gluon vertex by ([T]) and ([2]) of Fig. [TJ we need to introduce a truncation 
at some point in the infinite tower. This need for a truncation also applies 
to the BSE since therein we must specify the (2PI) four-point scattering 
kernel. This is a delicate process since we must be careful to preserve 
various symmetries of the theory. The most important of these relates to 
chiral symmetry, expressed via the axial-vector Ward-Takahashi identity 
(axWTI). This underpins the observed mass spectrum of the light mesons, 
and ensures that pions are indeed the (pseudo)-Goldstone bosons of the 
theory. This identification is a necessary feature of any serious model and 
must be exhibited by any truncation of the BSEs and DSEs. 

The simplest truncation that satisfies this criterion is that of Rainbow- 
Ladder (RL) whereby the full quark-gluon vertex is replaced by a bare 
vertex, see e.g. [1]. The axWTI preserving kernel in the BSE then corre- 
sponds to a single gluon exchange, re-summed to all orders thus providing 
the 'ladder'. These RL models are designed to reproduce predominantly 
s-wave mesons due to the simple vector- vector structure of the interaction. 
To compensate for this simplicity one constructs a phenomenological model 
for the gluon, effectively subsuming additional vertex corrections from the 
Yang-Mills sector and unquenching effects. 




Fig. 1. DSEs for: (p} fully dressed quark propagator; ([2]) full quark-gluon vertex; 
([3]) truncated quark-gluon vertex. Internal propagators are dressed, with gluons 
shown by wiggly lines, quarks by straight lines and dashed lines mesons. White- 
filled circles show meson amplitudes whilst black-filled represent vertex dressings. 
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To separate the phenomenological from the ab initio, we must investigate 
the quark-gluon vertex and determine the impact of corrections beyond tree- 
level to our quarks and mesons. Such studies have been made e.g. [2-8]. 
Following the analysis of [9, 10] we approximate the full DSE (|2|) with the 
(nonperturbative) one- loop structure of ([3]). Here the first 'non-Abelian' 
loop-diagram in ([3]) subsumes the first two diagrams in the full DSE to first 
order in a skeleton expansion of the four-point functions. We neglect the 
two- loop diagram in the full DSE ((2J), which is justified for small and large 
momenta [6,9]. The remaining 'Abelian' contributions are split into the 
non-resonant second loop-diagram in ([3]) and a third diagram containing 
effects due to hadron back-reactions. 

In the next section we consider these resonance contributions, using a 
RL truncation for the non-resonant parts. We follow this by exploring a 
new truncation in which leading non-resonant parts from the non-Abelian 
vertex are included self-consistently in the quark DSE and meson BSE. 

3. Including unquenching effects 

The prescription for including pion degrees of freedom in the DSEs and 
BSEs in a manner consistent with the axWTI have been proposed and 
investigated in [10-12], with the resultant system of equations depicted by 




The first loop diagram of both equations relates to the usual rainbow-ladder, 
where the infrared suppressed gluon is enhanced by the vertex dressing, 
indicated by YM. The second diagram contributing to the quark DSE and 
meson BSE represents the back-reaction of the pion onto the quark. This 
requires input of the quark-pion vertex, which we parameterise by a chiral 
approximation of the leading pion Bethe-Salpeter amplitude [12] 

T{(p;P)=r^ 5 B x (p 2 )/f n . (4) 

Here B x (p 2 ) is the scalar dressing function of the quark propagator in the 
chiral limit. The effects of neglecting the three sub-leading amplitudes have 
been quantified for a real- value approximation in Ref. [10], and found to be 
on the level of a few percent. The advantage of the approximation in 
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is that we can then directly calculate the quark propagator in the complex 
plane. 

We need to specify the gluon propagator and a quark-gluon vertex that 
subsumes the non-resonant parts of the interaction into an effective rainbow- 
ladder model. Since here we wish to employ a gluon propagator as calculated 
from its DSE, we employ the soft divergent model (SD) for the quark-gluon 
vertex as described in [9,12]. 

We calculated a range of meson observables and observe that the effect 
of the pion back-reaction has only a small impact on the pion mass itself. 
The impact of including pion-cloud effects on the leptonic decay constant is 
more substantial, with effects of the order of 10 MeV: 
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For the remaining heavier mesons, the common trend is that the inclu- 
sion of the pion cloud decreases the masses by 100-200 MeV. Most notable 
of these are for the rho, where we see that unquenching from the pion-cloud 
yields a bound-state that is ~ 90 MeV lighter. This will prove important in 
what follows when we consider gluon self-interaction contributions. 

It is clear, however, that in order to reproduce the rich spectrum of light 
mesons we need to include spin dependent contributions from the Yang- 
Mills part of the quark-gluon vertex. We now consider this in the absence 
of pion-cloud effects in the next section. 



4. Including gluon self-interactions 

Looking back to our proposed truncation of the quark-gluon vertex, we 
identified the second loop diagram of @ as the dominant contribution. Ig- 
noring all other contributions the resulting equation, coupled with the quark 
DSE, may be solved numerically provided we know the gluon propagator 
and three-gluon vertex dressing. Since we work in Euclidean space we need 
the quark propagator for complex values of its momenta. This means our 
quark-gluon vertex must also be solved for complex momenta. Through 
judicious choice of momenta in both the quark and vertex DSE, this can be 
performed without unconstrained analytic continuation of the gluon prop- 
agator and three-gluon vertex. The axWTI preserving truncation for the 
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BSE [2, 13], consistent with our choice of the quark-gluon vertex is [14] 




Since this work is a preliminary study of a new and sophisticated trun- 
cation [15] we choose to simplify things further by employing a simple mo- 
mentum dependent ansatz for the gluon [16] 

aZ(q 2 ) = (ttD/lo 2 ) q 4 e~<?l^ , (5) 

with two parameters D and u which provide for the scale and strength of 
the effective gluon interaction. Naturally, this ansatz provides only a first 
step towards a full calculation including input from the DSEs for the three- 
gluon vertex and the gluon propagator. Nevertheless we believe that the 
ansatz ([5]) is sufficient to provide for reliable qualitative results as concerns 
the effects due to the non-Abelian diagram onto meson properties. 

The results of our calculation follow, wherein we compare our truncation 
including gluon self-interaction effects in all twelve tensor structures of the 
quark-gluon vertex to that of RL with only vector- vector interactions. 
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The model parameters uj and D were tuned such that for the latter we 
obtain reasonable pion observables, and the quark mass is fixed at 5 MeV. 

What we find is that the mass of the p-meson is enhanced by ~ 120 MeV 
as compared to pure RL. This is intriguing since it has long been suspected 
that corrections beyond RL approximately cancel in the p [2,3,18]. This is 
supported by known estimates of mass shifts from the resonant and non- 
resonant Abelian diagrams in ([3]), calculated at ~ 90 MeV (see previous 
section) and 30 MeV [5] respectively 

5. Conclusions 

We presented a study of light mesons using two truncations of the Bethe- 
Salpeter equations beyond RL. These considered unquenching effects associ- 
ated with the pion cloud, and the dominant non-Abelian corrections to the 
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quark-gluon vertex stemming from gluon self-interactions. For the latter 
truncation, we obtain masses for the rho meson of ~ 900 MeV, consistent 
with quenched lattice simulations. The subsequent inclusion of unquenching 
effects and other non-resonant contributions to the quark-gluon interaction 
brings the rho mass back to its physical value thus supporting a long sus- 
pected cancellation mechanism. An investigation based on Yang-Mills DSE 
results together with unquenching effects is currently underway. 
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